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Qh| Abstract. The entanglement, purity and energy of two isolated two- level atoms 

which are initially prepared in Bell state and each interacts with a thermal cavity field 
are investigated by considering the atomic motion and the field-mode structure. We 
^ ' achieve the analytical solutions of the atomic qubits by using the algebraic dynamical 

O^. approach and the influences of the field-mode structure parameter, the strength of the 

thermal field and the detuning on the entanglement, purity and energy are discussed. 
We also investigate the state evolution of the atomic qubits based on the entanglement- 
JC . purity-energy diagrams. Our results show that the disentanglement process of the 

("^ I atomic qubits accompanies by excitations transferring from atomic subsystem to cavity 

VO ' field modes and atomic state from a pure state convert to the mixed states. 

rn 

in 
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^ ' 1. Introduction 

Entanglement is one of tlie most remarkable features of quantum meclianics and lias 
many practical applications in quantum information processing [Ij. However, realistic 
quantum systems are inevitably influenced by the surrounding environment, which 
always leads to decoherence of the quantum states. Particularly, the thermal field is 
frequently discussed in this problem. A thermal field, which is emitted by a source in 
thermal equilibrium at temperature T, is a highly chaotic field with minimal information 
about its mean value of the energy. However, such a chaotic field can entangle qubits that 
are prepared initially in a separable state [2], lead to entangled states in the interaction 
of a single qubit in a pure state with a thermal field regardless of the temperature of the 
field and reduce the system to a mixed state when the field variables are traced over [3] . 
The influence of the thermal field strength on the atom-atom entanglement |1] and 
atom-field entanglement [5] have also been investigated. Besides, Zheng |6] proposed 
a scheme for realizing two-qubit quantum phase gates with atoms in a thermal cavity. 
Jin [7] suggested a scheme of teleporting a two-atom entangled state with a thermal 
cavity and the success probability can reach 1.0. 



In this paper, we consider two isolated two-level atoms each interacting with a 
single-mode thermal cavity field. The effects of the atomic motion and the field-mode 
structure are considered at the same time. The atomic motion and the field-mode 
structure not only lead to nonlinear transient effects in the atomic population [SI [9], 
which are similar to self-induced transparency and adiabatic effects, but also give rise 
to the periodic evolution of the entropy squeezing [10], the field entropy, the atomic 
inversion [11] and the entanglement [H |5]. Other effects by regulating the field- mode 
structure parameters have also been observed, e.g. decreasing the squeezing in the two- 
photon JC model [12], operating the entanglement and realizing the quantum gate 
operation [13]. The recent cavity quantum electrodynamics experiments which use 
an atomic beam passing along the axis of a rectangular or cylindrical cavity provide 
the feasibility of discussing the different field mode structures [Ml |T5]. Then, in the 
present work, we are interested in the effects of the atomic motion and the field-mode 
structure on the time evolution of the atomic state in thermal cavity field environment. 
We suppose the two atoms are initially prepared in one of the Bell states, which is a 
maximal entangled pure state. The previous studies of maximally entangled states 
usually use entanglement, purity and energy to characterized the set of two-qubit 
states [m [T71 [HI [19] . As the atom- field interacting process accompanied by exchanging 
excitations between atoms and fields, the energy in the system is a direct infiuences on 
the entanglement and mixedness properties. On the other hand, the initial pure state of 
the atomic qubits must be infiuenced by interacting with a thermal cavity environment, 
so the mixedness is an important characteristic in the process of entanglement evolution. 
As a result, we will investigate the time evolution of the entanglement, purity and energy 
via manipulating of the field-mode structure parameters, the strength of the thermal 
field and the detunings between atoms and thermal cavity fields. The relationships 
among entanglement, pure and energy will also be presented with entanglement-pure- 
energy (EPE) diagram, which can offer a nice visual to the allowed state of the atomic 
qubits. 

The paper is organized as follows: In Sec. 2, we first describe the model under 
consideration and then derive the exact expression for the atomic reduced density matrix 
using algebraic dynamical approach [20l [211 1221 ES] . The quantities used to quantify the 
entanglement, purity and energy of atomic qubits are also defined in this section. Sec. 3 
is devoted to investigate the time evolution of entanglement, purity and energy for the 
atomic qubits. In Sec. 4, we discuss the time evolution of the atomic qubits with a EPE 
diagram. Finally, we present our conclusion in Sec. 5. 

2. Model 

We consider two identical moving two-level atoms (A and B) and two spatially separated 
cavities (a and b) with non-decaying single mode fields by using very high quality factor 
cavities [HI [2^ . Atoms A and B fiy through cavities a and b with a constant velocity, 
respectively. We suppose the two subsystems Aa and Bb are identical with same value 



of atom-field coupling strength, frequencies and field-mode structure. The Hamiltonian 
for the considered system in the rotating- wave approximation can be written as (/i = 1) 

H =Hi + H2, 

Hi = u,a^a + uoSf + gf{z){a^S^ + aS^), 

H2 = co^b^b + cuoSf + gf{z){b^S^ + bS^) (1) 

where Hi and H2 are the Hamiltonians for subsystems Aa and Bb, respectively, a^ 
and a (6"'" and b ) are the creation and annihilation operators of the cavity field a (b). 
S]_, 5*1 and SI represent the atomic raising, lowering and inversion operators of the 
atom i{i = A,B). Uc and uq are the frequencies for the field a and the atom A (or the 
field b and the atom B), respectively, g is the atom-field coupling strength, and f{z) is 
the shape function of the cavity field mode. When the interaction energy of atom-field 
coupling is much larger than the transverse kinetic energy spread of the atom, we can 
neglect the transverse velocity spread and restrict our investigation to atomic motion 
along the cavity axis {z axis). Then the atomic motion is incorporated into f{z) as 
follows 

/W ^ fivt), (2) 

where v is the atomic motion velocity. In this regard the cavity field-mode TEM^np is 
defined like f{vt) = sin{p7rvt/L), where p represents the number of half wavelengths of 
the field-mode inside a cavity with length L. If the atom passes through the cavity so 
fast that the atomic motion can be considered as a constant. For a proper choice of the 
atomic motion velocity v = gL/n, then j^ f{vt )dt = [1 — cos{pgt)]/pg. 

In the following, we propose the algebraic dynamical approach to derive the time 
evolution operator and the density operator based on the Hamiltonian (1). The key idea 
of the algebraic dynamical approach is introducing a canonical transformation operator 
that transforms the Hamiltonian into a liner function in terms of a set of Lie algebraic 
generators. According to algebraic dynamics, linear systems are integrable and solvable, 
then the time evolution operator and the density operator can be obtained easily. In the 
case of symmetric atom- field interaction, the two subsystems are completely equivalent. 
For simplicity, we will work with the subsystem Aa. A straightforward analysis of the 
Hamiltonian (1) shows that the total excitation number for subsystem Aa is 

Ni = a^a + Sf + ^, (3) 

which is a conserved quantity for the subsystem Aa and commutes with the Hamiltonian 
Hi. Based on the algebraic dynamical approach, introducing SU(2) algebra generators 
{Jo, J+, J-}, with Jo = Sf, J^ = N^ ' aS^, J_ = A''^' ' a^S^, which are nonlinear 
expressions and obey the following commutation relations 

[Jo, J+] = J+, [Jo, J-] = -J-, [J+, J-] = 2Jo. (4) 

In terms of the SU(2) algebra generators and the canonical transformation operator 
Ug = exp(6' J+ — ^ J_), we can obtained the time evolution operator of the subsystem Aa 

Ui{t) = e-*^^-^* = ?7,e'*(^«"^i^«)*[/;i 



= e-*-^i*[cos 2iJo sin — cos2^ + ?(J+ + J_) sin — sin 2^]. (5) 

where Ei = w,(iVi - i), ^ = - arctan[(v/AV4T7^ - ^/^)/9'nI^^], A = 
a/A2 + Ag''^Ni, g' = ga/t, a = f^ f{vt')dt' = [1 - cos{pgt)]/pg and A = uq - Uc is 
the detuning between the atom A and the cavity a. 

What should be noticed here is that using canonical transformation operator to 
diagonalize the nonlinear Hamiltonian (1) doesn't change its intrinsic qualities. Likewise, 
we can get the evolution operator U2 of the subsystem Bb, which has the similar form 
as Ui. 

f/2 (t) = e~^^2t [cos !_ _ 2iLo sin — cos 20 + i (L+ + L_ ) sin — sin 20] . (6) 

where E2 = Uc{N2 - i), = - arctan[(v/AV4T7^ - ^I'^)I9'nI% and r] = 
a/A^ + Ag'^N2. N2 = b'^b + S^ + | is the total excitation number for subsystem 
Bb. {Lq, L_|_, L_} are the SU(2) algebra generators with Lq = S^, L^ = N2 bS^, 

Throughout this paper we suppose the two atoms AB to be initially prepared in 
one of the Bell states, |\E') = -y^{\eg) + \ge)), and the two thermal cavity fields ab are in 

single-mode thermal field states Pa(0) = X]^o-^«l^)(^l' Pb{^) ~ X]m=o-^»nl"^)("^l- ^^ ^ 
result, the initial density operators for the two atoms and the two thermal cavity fields 
can be written as 

Pab{0) = |^)(^| = -{\eg){eg\ + \eg){ge\ + \ge){eg\ + \ge){ge\), 

00 00 
P/(0) =Pa(0)®p6(0) = ^^P„P„|nm)(nm|, (7) 

n=0 n=0 

where P„ = j^^^,Pm = ^i+^^+i ■ k = l/[exp(wc/Ta)-l] and / = l/[exp(a;c/T6)-l], k 
and / are the mean photon numbers of the thermal cavity field mode a and the thermal 
cavity field mode b, corresponding to the temperatures Ta and T^, respectively. 
Then, the initial density operator for the total system can be derived as 

pAB-f{0)=pAB{0)^PfiO) 

= — 2^ / ^ PnPmJlGiT'gm) {engm\ + \engm) {gnem\ 

n m 

+ Ignem) {engm\ + \gnem) (gneml) , (8) 

where the lengm) indicates that atom A is in the excited state and atom B is in the 

ground state, field mode A and field mode B are in the states \n) and \m), respectively. 

The initial state (8) under the action of the operator Ui{t) ® U2{t) evolves to 

PAB-f{t) = Ui{t)U2{t)pAB-f{0)Ul{t)Ul{t). (9) 

Then, from Eq.(9), we can get the reduced density matrix pABit) of the subsystem 
AB by tracing over the thermal cavity field variables. In terms of the atomic basis states 



\gg)^ Ifi'c), \^g), and |ee), the reduced density operator pAsif) can be expressed as 

/ xi \ 



pAB{t) = Trf[pAB^f{t)] 



a;2 Xs 
0:4 X5 

\ xe y 



(10) 



where Xi + X2 + x^ + Xq = 1, 
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and A„ = a/A^ + Ag''^n, rjm = a/A^ + Ag'^m, 9^ = - arctan[(A/A2/4 + g'^n - 
A/2)/g'n^/^], 0„ = - arctan[(v/AV4 + g'^m - A/2)/g'm^'\ 

Based on the analytical solution of pAB{t), we can conveniently do some 
approximation in analyzing the numerical results in the following section. Besides, 
we can define the quantities to quantify the entanglement, purity and energy of the 
atomic qubits and give their expressions in terms of the matrix elements of pAsit)- 
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We adopt Wootters' concurrence as a measure of entanglement in this 
discussion |25], which is denoted as Cab = Maa;{0, a/Ai — V^ — -y/As — V^}, and A, 
are the eigenvalues of the matrix (pabPab) in non increasing order. Following Eq.(lO), 
the expression of Cab turns out to be 



Cab = 2Max{0, \x3\ — a/xi x x^}. (12) 

When the value of Cab is positive, the atomic system is entangled. Cab = 1 
corresponds to the maximal entanglement state, while Cab = indicates the atom A 
and the atom B are separable. 

The energy Uab of subsystem AB is defined here as the expectation value of the 
Hamiltonian Hab = (^oS^ + ujqS^. We set coq = 1, then the energy Uab can be obtained 
based on the expression of Eq.(lO) 

Uab = Tr{pAB-f{t)HAB] = Xg - Xi (13) 

For two two-level atoms, Uab ranges from -1 for pab = \gg){gg\ to 1 for pab = 
\ee){ee\. 

To quantify the mixedness of the state PAB(t), we use the purity 

Pab = Triplsit)} = xj + xj + xj + xl + 2x3X4 (14) 

For the atomic qubits, Pab ranges from 1/d for completely mixed state to 1 for 
pure state for d-dimensional systems, which is closely to the linear entropy measure of 
mixedness. 

The aim of this paper is to address the question how the atomic motion and 
the field- mode structure infiuence the state pABif) of the atomic qubits in the cases 
of thermal environment. We know that the two atoms are initially in the maximal 
entanglement state, and the atomic state will evolves with time followed by the variation 
of the entanglement, purity and the transfer of the energy. Their time evolution will be 
discussed in the next section. 

3. Entanglement, purity and energy versus time 

There are three controllable parameters in the analytical expression of pAB{t)'- the field- 
mode structure parameter, the mean photon number in each cavity and the detuning. 
In this section, we will discuss their effects of the three parameters on the time evolution 
of the entanglement, purity and energy of the atomic subsystem. 

In Fig.l and Fig. 2 we plot the time evolution of Cab, Pab and Uab affected by 
different values of the field-mode structure parameters and the mean photon number in 
the situation of exact resonance. Form Eq.(ll), we can easily find that for resonant 
atom-field coupling, A=0, A„ = 2g'y/n, rjm = 2g'y/rn, sin26'„ = sin20m = —1, 
cos26n = cos 20m = 0. The elements of the matrix pab expressed in Eq.(ll) convert to 



Xi = -'^'^[Pn-lPmSm'^{g't^/n) cos^ {g't^/m) 



,.1 (al ,.1 <P) 




Figure 1. The evolution of Cab, Pab and Uab for different field-mode structure 
parameters, p — f(solid) and p — 4(dotted). The two subsystems are symmetric for 
k = I = Q.\ and the atom and the field are in exact resonance, (a) Cab versus time 
gt] (b) Pab versus time gt; (c) Uab versus time gt. 



+ PnPm~iCOS^{g'tVn)s\\i^{g't^/m)], 
X2 = - ^"^[PnPm cos^ {g't^/n) cos^{g'ty/m + l] 

n m 



+ Pn-iPm+1 sin^ig'ty/n) sin^{g'tV'm + 1)], 
X3 = - y^ y^[PnPm cos{g't^/n) cos{g't\/n + 1) cos{g'ty/m) cos{g't\/m + 1)], 

n m 

X4 = X3, 

X5 = - ^ ^[Pn+i-Pm-i sin^ {g'Wn + 1) siii^ig'ty^) 

n m 



+ PnPm cos^ {g'ty/n + 1) cos^ {g't^/m)], 
xe = - J2J2^PnPm+i cos^ig'Wn + 1) sin'^ {g'tVm + 1) 



2 

n m 



+ Pn+iPmSin^ig'Wn + 1) cos^ig'Wm + 1)]. (15) 



5". i. i/ie effects of field-mode structure parameters 

The atomic initial state is a Bell state, which is a pure state with P^^ = 1 as well 
as a maximal entanglement state with Cab = 1- However, the initial energy for the 
atomic subsystem is zero, that is Uab = 0. Fig.l illustrates the cases when the atom 
is in motion at the velocity v = gL/ir for parameters p = 1 and p = 4, respectively. 
It has been studied that when the atomic motion is considered, the time behaviors 
of field entropy, atomic inversion [11], and entropy squeezing [10] are periodical, and 
their evolution periods are shorten with the increase of parameter p. Similar behaviors 
occur in this work. From Fig.l we can find that the evolution periods is decreased 
with the increase of parameter p. This is because the time factor is the scaled time 
gt when the atomic motion is neglected, and is g't when the atomic motion is taken 
into account, g't = [1 — cos{pgt)]/p is a periodical function on the scaled time gt with 
period 2tt/p. In addition, the amplitudes for Cab, Pab and Uab are reduced while their 






Figure 2. The evolution of Cab, Pab and Uab for different thermal field strengths, 
k ^ I ^ 0.5 (solid), k = I — 5 (dotted). The atom and the field are in exact resonance 
and the field- mode structure parameter p — I. (a) Cab versus time gt; (b) Pab versus 
time gt; (c) Uab versus time gt. 



maximuni values are still unchanged. That is, compared with p = 1, the "sudden death 
of entanglement" disappears, the maximum mixedness of the atomic state reduces and 
the energy exchange between atoms and field modes decreases with the increase of the 
field-mode structure parameter p. 

3.2. the effects of mean photon number in each cavity 

What we talked about in Fig.l is just limited to the situation of very weak thermal 
field with mean photon number A; = / = 0.1 in each cavities. Then we are interested in 
how the atomic qubits evolves as the mean photon number in each cavity increases. It 
has been demonstrated that thermal cavity field can lead to entangled state of quantum 
qubits interacting with it |2l [3] , while strong thermal cavity field can inhibit the atom- 
atom entanglement |1] and atom-field entanglement [5]. Here, we pay our attention to 
the influence of mean photon number on the non-local atom-atom entanglement, in the 
case of exact resonance, as is shown in Fig. 2(a). In addition, in Fig. 2(b) and Fig.2(c), 
further study is employed to the time evolution of purity and energy, which can help 
us to get more information about entanglement evolution. From Fig. 2 we can find that 
with the increase of mean photon number in each cavity, both the amplitudes and the 
maximum values for Cab, Pab and Uab decrease. That is, compared with the case of 
weak thermal cavity field, the atomic qubits, which couples to strong thermal cavity 
fields, can not evolve to a maximal entanglement state with Cab = 1 as well as a pure 
state with Pab = 1- Meanwhile, the time interval of the "sudden death of entanglement" 
lengthens, the maximum mixedness increases, while the energy transfer between atomic 
qubits and field modes in each interacting period is more and more less. 



3.3. the effects of detunings between atom and field 

The effects of detuning on the entanglement, purity and energy between atom A and 
atom B are depicted in Fig. 3. In Fig. 3(a), Cab oscillates at first and as time evolves it 
will maintain the maximal entanglement state when the influence of the atomic motion 
is considered and A 7^ 0. Moreover, the larger the value of detuning, the faster the 





Figure 3. The evolution oiCAB, Pab and Uab for different detuning, A = O.l(solid), 
A ~ 1 (dotted) and A = 5 (dashed). The two subsystems are symmetric for /c = ^ = 0.1 
and the field- mode structure parameter p ~ \. (a) Cab versus time gt] (b) Pab versus 
time gt] (c) Uab versus time gt. 



Cab reaches the stable maximal value 1. When A ^ g, the oscillating almost vanishes 
and the two atoms nearly entangle maximally all the time as time evolves. Similar 
behavior to the time evolution of Pab is depicted in Fig.3(b). In Fig. 3(c) we can find 
that with the increase of A, the amplitude of Uab is reduced gradually, and large values 
of A make it more easy to get to zero. In a word, with the increase of detuning, 
the atomic subsystem is almost "frozen" in the initial state. This can be explained 
as follows: on the one hand, based on the expression of Pab, for weak thermal cavity 
fields A; = / = 0.1 and large detuning A ^ (7, A„ = 77™. ~ A, sin26'„ = sin20m ^ 0, 
cos 26n = cos 20„ ^1. As a result, a;i = xg ^ 0, X2 = X3 = X4 = ^5 ^ | Xln Zlm PnPm- 
Then, Cab = 2Max{0, \xs\ - V^i x x^} ^ E„E™^nP™ = 1, Pab = Tr{p%{t)} = 

Xl+ Xl + Xl + xl + 2x^Xi K. J2nJ2mPnPm = 1, UaB = Tr{pAB-f(t)HAB} = Xq-Xi ^ 0; 

on the other hand, the interacting process between atoms and fields is accompanied 
by the transfer of the excitation between the localized atom and cavity mode, which 
depends on the atom-field coupling and is distinctly infiuenced by the value of detuning. 
Larger detuning can inhibit the atom-field coupling and restrain this transfer process 
greatly, therefore the initially maximal entanglement pure state can be "frozen" in the 
atomic subsystem. 

4. Entanglement-purity-energy diagram 

In this section, we devote to investigate the relationships among entanglement, purity 
and energy for atomic subsystem, which refiects much of the nontrivial information 
about the particular atomic state in the atom-field interacting process. Here, we limit 
our study to the weak thermal field with k = I = 0.1 and the field- mode structure 
parameter p = 1 in resonant situation. Under these conditions, cos2^„ = cos20m = 0, 
sin2^„ = sin 20m = ~1- We plot entanglement-purity-energy diagram in Fig. 4(a) 
and show its projections on entanglement-energy and entanglement-purity planes in 
Fig. 4(b) and Fig. 4(c), respectively. At initial time. Cab = 1, Pab = 1 and Uab = 0, 
the atomic qubits in the maximal entangled state. From Fig. 4(b) we can find that the 
disentanglement process accompanies by excitations transferring from atomic subsystem 
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Figure 4. The relationship between entanglement, mixcdness and energy for p = 1, 
A ^ 0, k ^ I — 0.1. (a)three-dinicnsional diagram for Cab, Pab and Uab] (b) Cab 
versus Uab] (c) Pab versus Jab- 



to cavity field modes and atomic state from a pure state convert to mixed states. Tlie 
minimal energy for atomic subsystem is about -0.7 when the two atoms are separable, 
and the maximum value is zero when the two atoms are in the maximal entanglement 
state. This suggests the atomic state can not evolve to pab = \gg){gg\ {Uab = —1) 
or Pab = |ee)(ee| {Uab = 1) in the atom-field interaction process, which can also be 
confirmed in Fig. 4(c). When Cab = 0, Pab 7^ 1, the atomic qubits is in mixed state 
when they are separable. While Cab = 1 corresponding to Pab = 1, this indicates 
that the atomic qubits can just realize the maximal entanglement pure state but the 
maximal entanglement mixed state can not be obtained. 

5. Conclusion 

In this paper, we employed three parameters — entanglement, purity and energy to 
describe the information about two distant atoms which are initially prepared in Bell 
state. Our results showed that considering the atomic motion and the field-mode 
structure can lead to the periodic evolution of entanglement, purity and energy. With 
the increase of field-mode structure parameter p, both their evolution periods and their 
amplitudes are decreased while their maximum values are unchanged. However, strong 
thermal field can reduce the peak values of entanglement, purity and energy of the atomic 
qubits and make the atomic state initially in a pure state to mixed states. Meanwhile, 
in such a chaotic field, energy transfer between atoms and fields is more and more less 
with the increase of thermal field strength. In addition, large detuning is in favour of 
reducing their oscillation time and "frozing" the initial maximal entanglement state in 
the atomic subsystem. We also analyzed the possible state that the atomic qubits may 
evolves into. From the entanglement-purity-entanglement diagram we found that the 
disentanglement process for the atomic subsystem accompanies both by the excitation 
transferring from atomic subsystem to cavity field modes and the state converts from 
a pure state to mixed states. Our number results showed that, when the atomic state 
is in the maximal entanglement state, it is in a pure state at the same time; when the 
two atoms are separable, the atomic state is in a mixed state. However, in the atom- 
field interacting process, the state for atomic qubits can not evolve to the maximal 
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entanglement mixed state. 
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